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ABSTRACT 

We report the discovery of RG1M0150, a massive, recently quenched galaxy at z = 2.636 that is 
multiply imaged by the cluster MACSJ0150.3-1005. We derive a stellar mass of logM^ = 11.49lo‘i6 
and a half-light radius of i?e,maj = 1.8 ± 0.4 kpc. Taking advantage of the leasing magnification, 
we are able to spatially resolve a remarkably massive yet compact quiescent galaxy at ^ > 2 in 
ground-based near-infrared spectroscopic observations using Magellan/FIRE and Keck/MOSFIRE. 
We find no gradient in the strength of the Balmer absorption lines over 0.6i?e — 1.6i?e, which are 
consistent with an age of 760 Myr. Gas emission in [N II] broadly traces the spatial distribution 
of the stars and is coupled with weak Ha emission (log [N II] /Ha = 0.6 ± 0.2), indicating that 
OB stars are not the primary ionizing source. The velocity dispersion within the effective radius 
is cTe, stars = 271 ± 41 km s“^. We detect rotation in the stellar absorption lines for the first time 
beyond z ^ 1. Using a two-integral Jeans model that accounts for observational effects, we measure 
a dynamical mass of logMdyn = 11.24±0.14 and U/cr = 0.70±0.21. This is a high degree of rotation 
considering the modest observed ellipticity of 0.12 ± 0.08, but it is consistent with predictions from 
dissipational merger simulations that produce compact remnants. The mass of RG1M0150 implies 
that it is likely to become a slowly rotating elliptical. If it is typical, this suggests that the progenitors 
of massive ellipticals retain significant net angular momentum after quenching which later declines, 
perhaps through accretion of satellites. 

Subject headings: galaxies: elliptical and lenticular, cD—galaxies: kinematics and dynamics— 
gravitational lensing: strong 


1. INTRODUCTION 

Observations are beginning to reveal the bulk proper¬ 
ties of the high-redshift progenitors of massive elliptical 
galaxies. Massive quiescent galaxies are present at 2 : > 2 
but have very compact sizes compared to local early- 
type galaxies (ETGs). Much effort in recent years has 
focused on understanding the growth of these “nuggets” 
and the emergence of extended stellar envelopes (e.g., 
Trujillo et al.||2QQ6| Ivan Dokknm et al.||2QlQl [Newman 


111 _ 

et al.||2012| |Patel et al.||2013| |Belli et al.]|2015p'. Recent 

observations are also beginning to link these early qui¬ 
escent galaxies to star-forming progenitors at eve n ear¬ 
lier epochs via their masses, sizes, and kinematics ( Barro 


Iier epocns via tneir masses, sizes, and kinematics ( Dar 

et al.l[2M4 Toft et ar]|2Q14 van Dokknm et aL]|2015| . 


JNumerical simulations have made increasingly detailed 


mergers (e.g.. 

Gox et al. 20061 Robertson et al. 

2006: 

Hopkins et al. 

2UU<) Zolotov et al. 2ul5|. I'orming a 


such gas-rich mergers are predicted to produce rotating 
remnants—pos sibly with a velocit y exceeding the bulk of 
today’s ETGs (jWuyts et al.|2010|). Different channels for 
the formation of compact galaxies can be tested with stel¬ 
lar population gradients ( Wellons et al.|[2M5 ) and kine¬ 
matic data. However, sucn comparisons require spatially 
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resolved observations, which are very challenging given 
the small angular sizes of typical 2 : > 2 quiescent galaxies 
(^ 07l — 072). Gravitational lensing provides a promis¬ 
ing route toward resolving these compact systems. The 
main challenge in locating lensed examples is the need 
for wide-area near-infrared (NIR) imaging to overcome 
their rarity and faintness at optical wavel engths. Thus 


far only three example h ave been published (Muzzin et al. 
2Q12[ Geier et al.||2Q13 ), and none has spatially resolved 


spectroscopic data. 

In this Letter we present the discovery of RG1M0150, a 
rare example of a massive, quiescent galaxy at 2 : = 2.636 
that is multiply imaged by the cluster MACSJ0150.3- 
1005 (R.A. 1^5Q^2 U.3, Deck -10°05'30", z = 0.365; 
Ebeling et al.||2001 ). RG1M0150 was located in a search 


tor bright sources with red NIR colors that are magnified 
by clusters in the Hubble Space Telescope (HST) archive. 
It is extremely bright {Ks^ab = 19.2) and well suited for 
detailed follow-up observations. We analyze NIR spectra 
collected with Magellan/FIRE and Keck/MOSFIRE and 
are able to spatially resolve stellar absorption lines with 
spectral resolution sufficient for kinematic work. We use 
these unique data to study the resolved stellar popula¬ 
tions and rotational support of RG1M0150. 

Throughout we adopt AB magnitudes and the cosmo¬ 
logical parameters (U^,U^,h) = (0.3, 0.7, 0.7). 
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Fig. 1.— Left: F814W/F110W/F140W composite image of MACSJ0150.3-1005 with an arcsinh stretch. Images 1, 2, and 3 of 
RG1M0150 are marked. Right: Zooms on the observed F140W images (panel b, linear stretch) are compared to the best-fitting 
Sersic model that has been lensed and PSF-convolved (panel c) and the residuals (panel d = b — c). Panel e (arcsinh stretch) shows 
the PSF-deconvolved images in the source plane, constructed by adding the unconvolved image plane model to the residuals (i.e., 
Data — Modelconvolved + Modelunconvolved) cind casting pixels to the source plane. 


2. OBSERVATIONS 

2.1. Imaging 

RG1M0150 was identified in a search of archival 
i/Sr/WFCS images of massive clusters. We targeted 
bright, magnified sources with red NIR colors indica¬ 
tive of a Balmer/4000 A break at z > 2. HST imag¬ 
ing of RG1M0150 covers WFG3/F140W, WFG3/F110W, 
AGS/F814W, and WFCP2/F606W (SNAP-12884 and 
SNAP-11103, P.L Ebeling; see Figure [^. We also ob¬ 
tained J, and K,^ band obser vations with Magel- 
lan/FourStar (Persson et al. 2013 ) which were reduced 
following procedures discussed by Kelson et al. (2014). 
Golors were measured in a 0'.'6 x 0. 75 rectangular aper¬ 
ture along Ima ge 1 matching the spectroscopic aperture 
used in Section |2.2[ Matching of the point spread func¬ 
tion (PSF) between the Fo urStar and HST image s was 
performed as described by Newman et al. (2012). To¬ 
tal fluxes were obtained by scaling to the F140W flux 
in larger 8" x 275 rectangle. We also incorporate 3.4/im 
and 4.6/im measurements from the AllWISE catalog]^ 
Zeropoint and PSE matching uncertainties of 0.03 mag 
(HST), 0.05 mag (EourStar), and 0.1 mag (WISE) were 
added in quadrature. 

2.2. Spectroscopy 


simultaneously. The MOSEIRE data cover the H band 
only but afford somewhat cleaner sky subtraction owing 
to the faster dithering cadence. A position angle (PA) of 
—7.5° was selected to place the slit along Image 1. 

The MOSEIRE exposures total 260 minutes taken over 
2014 November 26-27, consisting of 2 minute integrations 
in an AB dither pattern. The 077 slit provided a resolu¬ 
tion of (Jinst = 36 km s“^. Another slit was positioned on 
a star to measure the average seeing of 0767. Data were 
reduced using the MOSEIRE Data Reduction Pipeline. 

The EIRE observations were conducted over 2014 
November 1 and 3 and total 390 minutes, consisting of 
30 minute exposures using up-the-ramp sampling. The 
6" X 0775 slit provided a resolution of ainst = 30 km s“^. 
The average seeing of 0755 was determined by monitor¬ 
ing a star on the NIR slit viewing camera. The data 
were reduced using a modified version of the EIREHOSE 
pipeline that produces two-dimensional rectified spectra. 
Eor both data sets, telluric absorption correction and flux 
calibration wer e performed using AOV star observations 
and xtellcor ( Vacca et al.||2003 ). 


2.3. Lens Model 

We construct a lens model using the three images of 
RG1M0150 as constraints. Analytic models describe the 
lens mass distribution and the source light profile. Their 


age 1 of RG1M0150 with both MOSEIRE (McLean et al. 

tracing code described in iNewman et al. 

(2015), which 

2012 

) at the Keck I telescope and EIRE (^Simcoe et al. 

hts the pixel-level WF'G3/F140W data. 


dispersed echellette instrument that covers 0.8 — 2.45 /am 

^ http://wise2.ipac.caltech.edu/docs/release/allwise/ 


was modeled using an elliptical Sersic profile. To model 
the cluster lens, we use several dual pse udo-isothermal el- 
liptical (dPIE) mass distributions (see [Elfasdottir et al. 
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Fig. 2.— Left: Source plane F140W surface brightness profiles reconstructed from each image of RG1M0150 (solid colored lines). Dotted 
colored lines show the corresponding reconstructions from the PSF-convolved Sersic model fits, with residuals indicated in the lower panel. 
(At small radii, the reconstructions differ because the PSF affects each image differently.) Shading indicates the Icr uncertainty. The black 
dashed line shows the unconvolved Sersic model. Right: The compact size of RG1M0150 fblac k) is compared to 2.2 < ^ < 3 field galaxies 
drawn from the 3D-HST survey catalogs ( [Skelton et al.|2014| ) with structural parameters from [van der Wei et al. ]( [2014l >. 


2007). Each is described by two scale radii Vcut and Vcore 
and a normalization ao, in addition to its center, PA, 
and ellipticity. All parameters of the dPIE describing the 
cluster dark matter halo are free other than Vcut^ which 
is beyond t he strong lensing zo ne and can be fixed to 
1 Mpc (e.g., Richard et aH 2 C)lQ ). Additional dPIE com¬ 
ponents account for deflection by luminous cluster galax¬ 
ies and are generally modeled using the scaling relations 
(To = and Tcut = r* with prior s 

on ctq and as described by | Newman et al. (2013). 
However, the brightest cluster galaxy and the two galax- 
ies circled in Eigure[^ which are important perturbers, 
were allowed to depart from these scaling relations. We 
also included external shear to fine tune the image recon¬ 
struction, although this did not make a major difference 
to our final model. 

Eigure[^ demonstrates that our best-fitting model pro¬ 
duces three images that match the E140W data well. 
The model produces areal magnifications of fii = 3 . 9^0 7 , 


fi 2 = 2.41q 5 and /is = d.d^Q g for Images 1, 2 , and 3, re¬ 
spectively. Uncertainties were conservatively estimated 


by comparing to an independent 

set of Lenstool models 

(Kneib et al. 1993 

Julio et al. 

2007) constrained only 

by the positions, el 

lipticity, and fluxes of the images and 


varying the mass model components. 


3. STRUGTURE AND STELLAR POPULATION 
3.1. Stellar Structure 

Our best-fitting Sersic model indicates that 
RG1M0150 is a compact, nearly round, concen¬ 
trated galaxy with an effective radius Re,maj = 
0723 ± 0705 = 1.8 ± 0.4 kpc (semimajor axis), Sersic 
index n = 3.5 ± 0.9, magnitude mFi 4 ow = 21 . 8 lo '25 
ratio hja = 0.88±0.06, and PA = 2 ° ±14°. Uncertainties 
were estimated from the scatter in parameters obtained 
when fitting the three multiple images individually 


(mFi 4 ow a.lso includes the systematic magnification 
uncertainty). 

We reconstruct source plane surface brightness profiles 
for each image in Eigure The lensing magnification 
enables us to probe down to scales of R = 700 pc = 0.4i7e 
before PSE blurring becomes significant. Over the range 
R = 0.4Re ^ single Sersic profile remains adequate 

to fit the E140W light with systematic residuals of < 10% 
(lower panel). The rest-frame colors of RG1M0150 {U — 
V = 1.69±0.06, V — J = 0.73±0.14) satisfy the criterion 
often used to select quiescent galaxies. The right panel 
of Eigure demonstrates that RG1M0150 is a compact 
galaxy, lying slightly below the mass-size relation defined 
by other quiescent systems at similar redshift. 


3.2. Stellar Population 

RG1M0150 presents a red continuum shape with 
strong Balmer absorption dominating the rest-optical 
spectrum. In Eigure we combine the broadband pho¬ 
tometry with our EIRE spectrum to con strain the stellar 
popu lation using our pyspecfit c ode ([Ne wman et al 


2Q14| . We adopt so lar metall i city 


Bruzual ^ Chariot 


models and a Chabrier ( 2QQ3[ ) initial mass tunc 


tion. The EIRE spectrum is extracted in a rectangu¬ 
lar aperture extending ±076 from the center of Image 1, 
which best approximates ±Re in the source plane. 

We infer a remarkably high stellar mass of logM^ = 
ll.49lo.i6 5 scaled to the source plane Sersic model flux 
and including the magnification uncertainty. Given the 
very high indicated mass, we simulate the probability dis¬ 
tribution of the most massive 2 < z < 3 quiescent galaxy 
lensed by a cluster in the HST/WFC3 archive (80 clus¬ 
ters, assuming a typical cross-section a = 200 arcsec^ for 
ma gnification p > 3 sirn ilar to MAGSJ0150.3-1005) using 
thejMuzzin et al. (2013) mass functions. A lensed galaxy 
witn log > 11.49 is indeed quite rare {p = 0.005). The 
probability increases to p = 0.02 if the threshold is low- 
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Fig. 3.— Left: T he FI RE spectrum of Image 1, resampled to 300 km s“^ bins for clarity, is compared to the stellar population model 
described in Section |3.2| Green points show the broadband photometry, which is displayed on an expanded scale in the inset. Resampled 
two-dimensional spectra are shown above. Right: Zoom on the H(a+[N II] region with the stellar model (red) subtracted. Blue lines show 
a Gaussian decomposition with [N Il]A6585/[N Il]A6550 fixed to 3. The observed spectrum has been scaled by a factor of 3.6 to match the 
total fiux of Image 1. Gray bands denote strong sky lines. 


ered by Icr to logM^ > 11.33. If we take our simulated 
probability distribution as a prior in order to account for 
an Eddington-type bias, the posterior for the stellar mass 
would shift downward to logM>^ = 11.33lo!i5- 

For an exponentially declining star formation history, 
we find an age of 760 ± 50 Myr and a prompt e-folding 
time of 97 ± 11 Myr. The corresponding ongoing star 
formation rate is only SFR = 2 ± 1 Mq yr“^ (accounting 
for magnification), indicating that star formation was re¬ 
cently and quickly truncated in this galaxy. The specific 
star formation rate is logSSFR/yr“^ < —11.1 (68% CL) 
independent of magnification. Dust atten uation is mod¬ 
erate w ith Ay = 0.52 ± 0.08 assuming the |Calzetti et al.| 
(2000) reddening curve. (These uncertainties are purely 
statistical and do not include known differences among 
stellar population models.) 

The lensing magnification presents a unique opportu¬ 
nity to probe the location of the last star formation event 
by spatially resolving the Balmer absorption. We ex¬ 
tract spectra in 5 rectangular apertures whose bound¬ 
aries along the spectroscopic slit are chosen to exam¬ 
ine the same luminosity-weighted radius in the source 
plane on either side of the galaxy: (R) = 0.6Re: O.QRe, 
and 1.6Re- Figure demonstrates that the FIRE and 
MOSFIRE spectra agree well, and we therefore combine 
them in the bottom of the panel to search for gradients 
in Balmer line strength. No variation is apparent, even 
though the spectra would be sensitive to an age differ¬ 
ence of ^ 200 Myr. Thus, the last star formation event 
in RG1M0150 was not highly centrally concentrated. 

3.3. Ionized Gas 


We detect emission in the [N II] AA6550, 6585 and Ra 
lines (see Figure]^. Weak Ra emission with a modest 
equivalent width of EWr^ = 6 ± 2 A (rest frame) is evi¬ 
dent only after subtracting the underlying stellar absorp¬ 
tion. [N II] is much stronger (EW[nii]a6585 = 22 ± 2 A). 
The very high ratio log([N II]A6585/H(a) = 0.6 ± 0.2 ex¬ 
ceeds that attainable through photoionization by a star- 
burst (e.g., Kewley et al. ]20^ . The main ionizing source 
is therefore not OB stars, a conclusion consistent with the 
low star formation rate inferred from the stellar contin¬ 
uum. Although other bright emission line diagnostics are 
not available (e.g., [O II] and [O III] are inaccessible due 
to atmospheric opacity), similarly high [N II]/Ha ratios 
are sometimes seen in red sequence and post-starburst 
galaxies at lo wer redshift that usually have LINER -type 
spectra (e.g.,|Yan et al.||2006t |Lemaux et al.||2010|). The 
[N II]A6585 emission m RGIMUI5U has a roughly con¬ 
stant EW across the galaxy, which suggests (but does not 
prove) that the main ionizing source is distributed simi¬ 
larly to the stars (e.g., post-AGB stars or shocks) rather 
than an active galactic nucleus (AGN). X-ray, radio, and 
mid/far-IR observations could clarify this by testing for 
the presence of an AGN or obscured star formation. 

4. STELLAR KINEMATIGS: ROTATION IN A QUIESGENT 
Z >2 GALAXY 

Stellar kinematics were measured using the iL-band 
(^rest = 4060 A — 4930 A) sp ectra within the spatial 
apert ures described in Section |3.2[ We use the ppxf 
code (Cappellari & Emsellem||2004|) and allow a linear 

- TStf 


combination of Bruzual ^ Chariot 
populations as the template 


2003) simple stellar 
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4.1. Velocity Dispersion 
We first determine the velocity dispersion oy within 


Fig. 4.— Left: Resolved Balmer absorption lines extracted in several spatial bins. Corresponding spectra on opposite side of the galaxy 
center have been shifted to the systemic redshift and coadded. The FIRE and MOSFIRE spectra are combined in the lower spectrum 
and demonstrate no significant gradients in Balmer line strength over 0.6Re — 1.6Re- Right: The top panel illustrates the position of 
the spectrograph slit as seen in the WFC3/F140W image. The middle panel compares the fiux distributions in the spectra to the PSFs 
(dotted), demonstrating that we spatially resolve Image 1. Five extraction bins are indicated. The lower panel indicates the presence of 
rotation measured using the stellar absorption lines. 

oblate Jeans dynamical model. We assume mass fol¬ 
lows the light distribution described by our Sersic model. 
The free para meters are M^yn, the rotati on parameter k 
dSatoh 1980), and th e inclination. (See |van der Wei & 
van der Marel |2008 for a quite similar approach.) T'he 
dynamical model is constructed in the source plane, then 
lensed, convolved by the seeing, and binned like the data. 
The model is constrained by the rotation points (Fig¬ 
ure 13 and cTe, stars- Fortuitously, the slit PA of —7.5°, 
whi^ is aligned with the direction of magnification, is 
close to the major axis PA (2° ± 14°). 

The inferred dynamical mass is log M^yn = 11.24±0.14 
when marginalizing over the unknown inclination 0 < 
cosi < cosimin- Here we choose imin = 29° because it 
corresponds to a maximum intrinsic ellipticity eint =0.8. 
Although the intrinsic {V/a)[nt is sensitive to inclination, 
the projected value that we measure is not. The inferred 
Mdyn is also insensitive to inclination except when the 
galaxy is close to face-on (this result is gener ic in Jeans 
dynamical modeling, e.g 


the aperture encompassing (see Section 3.2). This 
is simply the second moment of the line-of-sigh^t veloc¬ 
ity distribution and includes both random and stream¬ 
ing motions. The FIRE and MOSFIRE spectra yield 
consistent estimates of (Je,stars = 288 ± 27 km s“^ 
and 254 ± 25 km s“^, respectively, which we average 
to obtain (Jg,stars = 271 ± 18 km s“^. Tests of sev¬ 
eral systematic uncertainties revealed that sensitivity to 
template construction and the inclusion or exclusion of 
Balmer lines dominate the error budget, contributing 
±38 km s“^. Combining the uncertainties in quadrature 
yields erg, stars = 271 ±42 km s“^, which is comparable to 
that measured for the ionized gas in the same aperture, 
= 255 ±32 km s-h 




e,gas 


4.2. Rotation and Dynamical Model 

Although our spectra are insufficient for resolved mea¬ 
surements of cr, we do have sufficient signal to measure 
velocities in our five spatial bins. This the first such mea¬ 
surement derived from stellar ab sorption features at red- 
shifts beyond the z ^ 1 study of [van der Wei fc van der 
Marel ( 2008|) . Eigure|^ (right panel) shows that all 4 in¬ 
dependent measurements on the southern half of the arc 
are redshifted, while all 4 measurements on the northern 
half are blueshifted with respect to ^^sys = 2.6356±0.0002 
measured in the central spatial bin. The presence of rota¬ 
tion with the same directional sense in independent data 
sets gives us significantly higher confidence in its reality. 

To compare the rotational velocity of RG1M0150 with 
other samples, we must account for observational effects, 
particularly the seeing, slit width, and binning. We 
do this by constructing a two-integral (semi-isotropic) 


van der Marel|1991 ). In partic¬ 
ular, when fixing i = imin we hnd that that the inferred 
Mdyn would increase significantly by 0.16 dex. Adopt¬ 
ing any inclination i > 32° that is even slightly larger, 
however, affects M^yn by < 0.05 dex. 

Adopting a local virial calibration log Mvir = 
log5crgRg/G = 11.19 ± 0.12 would have given a slightly 
smaller mass. M^yn is smaller than, but consistent with, 
our estimate of M^. However, some tension would arise 
had we adopted a Salpeter-type initial mass function. It 
will be interesting to evaluate this further with larger 
samples of lensed galaxies. 

Observational effects modify the observed rotation sig¬ 
nificantly, as Eigure demonstrates. Our modeling al¬ 
lows us to account for these. The best-fitting model has 
a maximum projected rotation velocity Wiax = 189 ± 
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Fig. 5.— Rotational suppo rt in RG1M0150 is compared to ob- 
servations of ETGs at 2 : ^ 1 (|van der Wei Sz van der Marel|2008| > 
and ^ Q (jEm sellem et al.|2011 ATLAS^^). The solid line is the 
|Binney| (|19Vy|) curve tor isotropic oblate rotators. The gray region 
appr oximately encompas ses the simulated high- 2 ; merger remnants 
from [Wuyts et al.| ( |2Q1Q| > . 

34 km and Via = 0.70 ± 0.21. This model adopts 
the classic definition of V/a based on Tmax and a^ aver¬ 
aged within i7e/2 (e.g., Davies et al.||1983 ). 


5. DISCUSSION 

A central question in the evolution of ETGs is the 
extent to which the quenching of star formation and 
the transformation of stellar structure, morphology, and 
kinematics proceed together (e.g., in one major merger 
episode) or in several distinct events. The presence of 
compact quiescent galaxies at high redshift demonstrates 
that growth of the stellar distribution proceeds long af¬ 
ter quenching. Recent examinations of the evolving el- 
lipticities of quiescent galaxies suggest tha t these early, 
compact systems were also mo re disk-like (van der Wei 
et al.||20H |Chang et al.||2013|). The key test is a direct 
measurement of rotation, presented in this Letter for the 
first time. 

Figure compares the rotational support of 
RG1M0150 with that of local ETGs in the classic 
(e, E/cr) diagram]^ RG1M0150 has a high V/a for its 
modest observed ^lipticity that places it on the upper 
envelope of local ETGs. Although we cannot determine 
the inclination and deproject V/a with current data, 
it is possible that the edge-on rotation velocity is 
significantly higher. 

Numerical simulations of galaxy mergers require high 
gas fractions of ^ 40% (Wuyts et al. 2010) to produce 
compact remnants such as RGllVlOloCr Tnese simula¬ 
tions generally find that a significant amount of angular 
momentum is retained and suggest that rotation will be 


more significant in quiescent galaxie s at higher redshift , 
when mergers were more dissipative. Wuyts et al. (2010) 
predicted that quiescent z ^ 2 galaxies would tall on the 
upper envelope of local ETGs in the (e, E/a) diagram. 
Intriguingly, we find that RG1M0150 indeed lies close in 
Figure 5 to their simulated merger remnants, which are 
compact systems with M* 10^^ Mq and ^ 1 kpc 
comparable to ty pical z ^ 2 quiescent galaxies. 

Simulations by Zolotov et al. (2015) predict that the 
gas in compact quiescent galaxies and their immediate 
star-forming progenitors is kinematically colder than the 
stars, i.e., has higher V/a. This causes the projected 
c^e,stars/c^e,gas to be ^ 1 for an cdgc-on view but pro- 
gressively hi g her fo r lower inclinations. Supporting this. 


Barro et al. (2015) measured cTgtars/c^gas = 1-7 ± 0.5 in 
a quenching compact galaxy at z = 1.7. We find that 
the stellar and gas line widths are closer in RG1M0150, 
with cTg stars/c^e,gas = IT ±0.2, but larger samples will be 
needed to test this hypothesis due to the strong inclina¬ 
tion dependence. 

Nearly all local ETGs with > 10^^-^ Mq (scaled 
to a Ghabrier IMF) are “slow rotators” with V/a much 
smaller than RG1M0150 (red points in Figure [^. Our 
best estimate for the stellar mass of RGlM015(Jalready 
places it within this range. Even allowing a somewhat 
smaller mass (e.g., our M^yn), it is quite likely that 
RG1M0150 will evolve into a slow rotator: empiri cal esti¬ 


mates based on number dens ity arguments (e.g., Muzzin 
et al.|2013 Patel et ^^[20 13 ) imply that RG1M0150 will 


grow m stellar mass by 0.2-0.4 dex by 2 : = 0. Stellar pop¬ 
ulations in such ultra-ma ssive galaxies today are almost 
exclusively very old (e.g., McDermid et al. |2015 ), so sig¬ 
nificant later star formation is unlikely, if kGlM0150 is 
typical, this suggests that the progenitors of giant ellip¬ 
ticals are still rotating significantly after they have been 
quenched and must reduce their net angular momentum 
later. A sequence of mostly dry mergers, suspected of 
driving the observed 4x growth in effective radius, 
might provide a physical explanation. 

General conclusions about galaxy formation obviously 
cannot be drawn from RG1M0150 alone. Fortunately, 
a larger sample of lensed compact galaxies at 2 : ^ 2 
(mostly at lower masses) is within reach and provides 
the most feasible route toward resolving the stellar pop¬ 
ulations and dynamics of these intriguing objects. 
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